Competition between final-state and pairing-gap effects in the radio-frequency spectra 

of ultracold Fermi atoms 
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The radio-frequency spectra of ultracold Fermi atoms are calculated by including final-state in- 
teractions affecting the excited level of the transition, and compared with the experimental data. 
A competition is revealed between pairing-gap efi^ects which tend to push the oscillator strength 
toward high frequencies away from threshold, and final-state effects which tend instead to pull the 
oscillator strength toward threshold. As a result of this competition, the position of the peak of 
the spectra cannot be simply related to the value of the pairing gap, whose extraction thus requires 
support from theoretical calculations. 

PACS numbers: 03.75.Ss,03.75.Hh,32.30.Bv,74.20.-2 
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One key issue concerning ultracold trapped Fermi 
atoms with a mutual attractive interaction is the detec- 
tion of their superfluid behavior below a critical temper- 
ature Tc- As conventional techniques (like the Meissner 
effect) to identify the occurrence of the superconducting 
phase do not work for these neutral systems, additional 
physical effects need be explored to the purpose. 

In this context, radio-frequency (RF) spectra were first 
proposed [2] to signal the presence of a pairing gap for 
ultracold Fermi (^Li) atoms at low enough temperature. 
Collecting data at various couplings and for temperatures 
from below to above Tc, the evolution of these spectra 
was qualitatively interpreted in the light of an underly- 
ing pairing gap. These experiments have since been ex- 
tended to different systems Q as well as in the presence 
of spin population imbalance [3|[, and quite recently also 
by attaining spatial resolution [J|. 

Several theoretical attempts have been put forward to 
interpret the RF spectra of ultracold Fermi atoms 0, 
SJl, S, i, K El- There is at present wide consensus 
[12| . however, that this difficult problem is still unsolved 
theoretically, and a challenging call has been launched 
[3| for a comprehensive theoretical interpretation of RF 
spectra of Fermi gases in the strong-interacting regime. 

A definitely important step forward for interpret- 
ing these spectra is the inclusion of final-state effects, 
whereby the atom excited to the final state |3) from the 
initial state |2) interacts with the atom in the state 
Here, the states |1) and |2) participate to the Cooper 
pairing as modulated by the Fano-Feshbach resonance 
with scattering length ai2, while the (empty) final state 
1 3) excited by the RF transition interacts with the state 
|1) via a different Fano-Feshbach resonance with scat- 
tering length ai3 0. This final-state interaction should 
considerably affect the density of states available to the 
excited atom in channel 1—3, thereby triggering a com- 
petition with pairing-gap effects present in channel 1—2. 

Purpose of the present paper is to show that consider- 
ing this competition considerably improves the descrip- 
tion of the RF spectra, by accounting for important fea- 



tures that cannot be explained otherwise. We limit here 
to low temperatures where the pairing gap is well devel- 
oped, but extension proves feasible also above Tc where 
the pseudogap takes the place of the pairing gap. 

On physical grounds, the approach discussed in the 
present paper bears strong analogies with the treatment 
of the exciton effect in the optical spectra of semiconduc- 
tors [l^. In that case, the long-range Coulomb attraction 
considerably modifies the absorption spectrum via the 
density of final states, leading to a finite (in contrast to a 
vanishing) value of the absorption coefficient at threshold 
[3]. For the RF spectra of ultracold Fermi atoms, final- 
state effects are less dramatic owing to the short-range 
nature of the interactions in both 1—2 and 1—3 chan- 
nels (which are assimilated by attractive delta-functions 
with scattering lengths ai2 and ai3, respectively). Yet, 
we shall see that for large enough jaial these exciton-like 
effects modify the RF spectra in a definitive fashion. 

The importance of taking the effect of 013 into ac- 
count was demonstrated in Ref. Q for the molecular 
(two-body) case, which represents the extreme BEC limit 
of the BCS-BEC crossover. In addition, Ref. con- 
sidered the effect of aia on the /-sum rule within the 
BCS-RPA approximation, albeit only in the BCS limit 
of the BCS-BEC crossover. In the present paper, we 
adopt a diagrammatic representation of the correlation 
function relevant to the calculation of the RF spectra 
of ultracold Fermi atoms, to include final-state on top 
of pairing-gap effects. This representation recovers the 
results of Refs. ^ and @ in the respective hmits. 

Besides the exciton-like effect between the photo- 
excited atom and the Cooper-pair-mate left behind, t- 
matrix self-energy corrections of level |3) due to interac- 
tion with level |1) may account for additional final-state 
effects associated with the action of the medium on the 
photo-excited atom. These self-energy corrections, how- 
ever, should not be included if atoms in level |3) do not 
reach thermal equilibrium with the environment in the 
time scale of the experiment. In addition, pairing fluc- 
tuations beyond mean field [3l are needed for a fuller 
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description of the initial state of the RF transition. 

Self-energy corrections of level |3) due to interaction 
with level |2) are instead dropped altogether, because 
they are irrelevant in the molecular limit and the inter- 
action 2—3 does not affect the value of the /-sum rule 
whenever level |3) is empty [§]. As the presence of a t- 
matrix self-energy in channel 2—3 would spoil this prop- 
erty, inclusion of additional diagrams for the response 
function would be required to recover it. 

In RF experiments with ^Li atoms, a rf pulse with 
frequency ul transfers population from state |2) to state 
|3) and the atom loss from state |2) is measured. For 
weak excitations, the RF signal 5{I{uj)) can be calculated 
within linear-response theory yielding [7]: 



j dvdv' Im{n^(r,r';( 



0} (1) 



where 7 contains the atomic matrix element of the tran- 
sition and n^(r,r';a;) is the Fourier transform of the 
retarded correlation function n^(r, r';i — t') = —iO{t — 
t'){[B{r,t),B{Y' ,t')\) taken at the frequency uj — ujl + 
/i — /13. Here, {r) = 'ip\{v)ip2{i^) is the relevant local 
transition operator and (• • •) a thermal average, while /i 
and /i3 are the chemical potentials for the initial (popu- 
lated) and final (empty) levels, in the order. 

Expression ([T]) applies to alternative dynamical ap- 
proximations. Following a standard approach for quan- 
tum many-body systems fl6| . the retarded correlation 
function is calculated via its Matsubara counterpart 

n(r,r';r!,) = / dr e*f^''^(T,[i^2(r', 0)Vl(r, t+) 

X V3(r,T)V'l(r',0+)]) (2) 

where T,- is the time-ordering operator for imaginary 
time T and = 2'kv/ (3 {v integer) is a bosonic Mat- 
subara frequency with inverse temperature (3. Analytic 
continuation from i^Li, to lo + i-q (where 77 is a positive 
infinitesimal) is taken eventually. 

Previous calculations 0, 0, Q attempting comparison 
with experiments considered a mean-field factorization of 
the thermal average in Eq.Q, resulting in the approxi- 
mate expression 



n(r,r';17,) 



- G2(r', r; u:^)G^{v, r'; + n,) (3) 



where aj„ = (2n -I- 1)7t/P {n integer) is a fermionic Mat- 
subara frequency and the G's are single-particle propa- 
gators for the initial |2) and final |3) levels of the tran- 
sition. This expression is represented by the "bubble" 
of Fig. 1(a). The RF spectrum calculated within the 
approximation fS]) shows the following features: (i) In 
the absence of interactions, the spectrum reduces to a 
Dirac delta-function 2tt^'^N25{ujl ~ ^a) where uia is the 
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FIG. 1; Diagrammatic representation of the correlation func- 
tion (O: (a) BCS bubble and (b) BCS-RPA series. Full and 
broken lines stand for fermionic single-particle propagators 
and interactions, respectively, while numbers attached to the 
end points identify the relevant states. The wiggly lines repre- 
sent the coupling to the RF field. For simplicity, all diagrams 
are drawn in four-momentum space. 



frequency of the atomic transition and N2 the total num- 
ber of atoms originally in state |2); (ii) Pairing interac- 
tions between states |1) and |2) shift this delta peak to 
higher frequencies and broaden it to the extent that it 
behaves like w"'^/^ for frequencies larger than the rele- 
vant fermionic energy scales, while preserving the sum 
rule for the total area of the spectrum [17 1 
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dLuS{I{Lu)) = 2Tr-f'^N2 



(4) 



(we assume throughout that level |3) is empty). 

The slow convergence for large uj resulting from the 
approximation ([3]) appears evident when one attempts 
comparison with the experimental spectra [10,11. It 
further leads to a divergent value for the /-sum rule that 
corresponds to the first moment of the spectrum i5(/(a;)). 
In this respect, it was pointed out in Ref. [§| that a fi- 
nite value for the /-sum rule results instead by going be- 
yond the bubble approximation ([3]) for i5(/(a;)) and sum- 
ming the BCS-RPA series, in analogy to the theory of 
superconductivity when im plem enting gauge invariance 
for the response functions [l8|. In Ref. [9|] this exten- 
sion was regarded appropriate to the weakly-interacting 
(BCS) regime of the BCS-BEC crossover, where compar- 
ison with the experimental spectra is hardly possible. 

Figure 1(b) shows the diagrammatic representation 
of the BCS-RPA series [l^ contributing to the corre- 
lation function in terms of which one can gain in- 
sights about the inclusion of the physical effects discussed 
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above, namely: (i) The occurrence of pairing in chan- 
nel 1-2 (the pairing-gap effect), which can be dealt with 
even by including fluctuations beyond BCS mean-field; 
(ii) The presence of an attractive interaction in channel 
1-3 (the exciton-like effect), which is evident in the lad- 
der Fis; (iii) The possible inclusion of the t-matrix self- 
energy dressing level |3) (the medium effect). [Pairing 
fluctuations and t-matrix effects need not be included in 
the internal lines defining Ti^.] Note that both pairing- 
gap and exciton-like effects are also borne out in the two- 
body calculation of Ref. @| . 

We have indeed verified that the results of the two- 
body calculation of Ref. Q can be recovered from the 
many-body calculation of S{I{lu)) based on the diagrams 
of Figs. 1(a) and 1(b), when approaching the BEC limit 
1 <C {kFCii2)~^ where kp is the Fermi wave vector asso- 
ciated with 7V2. By this remark, we can rely on the dia- 
grams of Figs. 1(a) and 1(b) to describe the RF spectrum 
d{I{uj)) throughout the whole BCS-BEC crossover, the 
only limitation being the condition of low-enough tem- 
peratures such that the pairing gap is well developed. 

The exciton-like effect contained in the diagrams of 
Fig. 1(b) generates a strong suppression of the tail of RF 
spectra, by augmenting its power-law convergence by one 
unity to a;~^/^. In addition, the diagrams of Fig. 1(b) 
give a vanishing contribution to the sum rule (jl]) (which, 
otherwise, could not be preserved). This implies that, to 
the negative tail for large lu resulting from the diagrams 
of Fig. 1(b), there corresponds a positive pile up for small 
LU which tends to peak the spectrum toward threshold. 

Figure 2 presents the results of our numerical cal- 
culation that includes pairing-gap (with fluctuations), 
exciton-like, and medium effects at the unitary limit 
(^^012)"^ = 0, and compares them with the tomographic 
spectra of Ref. Q. Good agreement between the theo- 
retical and experimental spectra results from this com- 
parison especially when the self-energy of level |3) is not 
included, thus suggesting that atoms in level |3) do not 
thermalize in the course of the experiment. In particular, 
our calculation is able to capture the main features of the 
experimental spectra, namely: (i) The value of the peak 
position normalized to the local Fermi energy, which in- 
creases only slightly from the trap center to the edge; (ii) 
Their asymmetry, the line shape being sharper on the 
low-frequency side of the peak; (iii) The strong progres- 
sive broadening of the spectra away from the trap center, 
which is accompanied by a lowering of the peak height. 

Note that the overall profiles of Fig. 2 reveal at glance 
the importance of the scattering length 013 for the RF 
spectra, since at unitarity for 012 spectra corresponding 
to different local densities present a non-universal behav- 
ior that can only originate from the presence of an addi- 
tional length scale in the problem. Note further that the 
first moment of these spectra (identified by an arrow in 
each panel) does not relate to the position of the experi- 
mental peak, being moved to higher frequencies owing to 
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FIG. 2: Comparison between theoretical [without (full lines) 
and with (broken lines) inclusion of the self-energy of level |3)] 
and experimental (dots) RF spectra for {kFai2)~^ = and 
(kpari)'^ = —0.93 (upper panel), —1.03 (middle panel), and 
— 1.32 (lower panel). These panels correspond to Figs. 2(b), 
2(c), and 2(d) of Ref. [ij, in the order, with radial position 
r — 0, 0.36, and 0.59 in units of the Thomas- Fermi radius, 
and associated temperatures 0.042, 0.052, and 0.085 in units 
of the local Fermi temperature. The value 013 = —174.6 nm 
is taken from Ref. [23[, while the local Fermi wave vector is 
taken from Ref. [j. The frequency is in units of the local 
Fermi energy and its zero is taken at the atomic transition. 
The spectral intensity is normalized so that the right-hand 
side of the sum-rule (4) is unity, while the arbitrary scale of 
the experimental data is fixed to match the theoretical peak 
height in the lower panel. 



the prominent role of the tail of the spectrum. 

To appreciate how the different physical effects con- 
tribute to the theoretical spectra of Fig. 2, we show 
in Fig. 3 how these spectra evolve with the progressive 
refinements of the approximations (for definiteness, we 
show the results for the lower panel of Fig. 2). We begin 
in Fig. 3(a) with the BCS bubble of Fig. 1(a), evaluated 
with the BCS gap and chemical potential. We further in- 
clude in Fig. 3(b) the BCS-RPA series of Fig. 1(b), with 
the same BCS values for the gap and chemical potential. 
We return in Fig. 3(c) to the bubble of Fig. 1(a), where 
we now include pairing fluctuations beyond mean field 
[lil in the fermionic propagator G2. Finally, in Fig. 3(d) 
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FIG. 3: Evolution of the theoretical spectra with the pro- 
gressive refinement of the approximations. In each plot, the 
experimental data have been normalized to match the theo- 
retical peak height. (See text for details.) 



we report the full calculation with all diagrams of Fig. 1 
included, and with pairing fluctuations consistently in- 
cluded in all fermionic propagators whose level |2) cou- 
ples to the external field (while no self-energy correction 
is considered for the fermionic propagator G3). In both 
Figs. 3(c) and 3(d) the gap and chemical potential are 
evaluated, too, with the inclusion of pairing fluctuations. 
In all cases, numerical accuracy results in the sum rule 
(HI) being satisfied at worse within 5%. 

The competition between pairing-gap and final-state 
effects on the RF spectra can be traced by examining 
in detail the above evolution. The simplest approxi- 
mation of the BCS bubble shown in Fig. 3(a) contains 
pairing-gap effects, which shift the peak of the spectrum 
to higher frequencies with respect to the atomic tran- 
sition, but does not include any final-state effect. The 
exciton-like effect is evident in Fig. 3(b), whereby the 
oscillator strength is pulled toward threshold with re- 
spect to Fig. 3(a) and the tail is strongly suppressed 
correspondingly. The role of the pairing gap is further 
emphasized by the inclusion of pairing fluctuations be- 
yond mean field in Fig. 3(c), which give rise to additional 
pseudogap effects that suppress the density of states near 
threshold and broaden the low-frequency side of the peak. 
The competing role of the exciton-like effect is again op- 
erative in Fig. 3(d), resulting in the RF spectrum that 
matches best the experimental data. 

We can finally compare the value of the peak position 
of the calculated RF spectra with the value of the pair- 
ing gap for a homogeneous system at unitarity. For this 
comparison we consider the full line in the upper panel of 
Fig. 2 at the lowest temperature, from which we extract 
the value 0.38 for the peak position (normalized to the 
local Fermi energy), while for the pairing gap Ref. 15 1 



provides the value 0.53 (the pairing gap extracted from 
the single-particle spectral function and order parame- 
ter about coincide for these coupling and temperature). 
This comparison confirms the role of the competition of 
the physical effects that we have discussed, leading us to 
conclude that, although RF spectroscopy undoubtedly 
provides a direct signature for the presence of a pairing 
gap in ultracold Fermi atoms, the numerical value of this 
quantity cannot be simply extracted from the peak posi- 
tion of the RF spectra. 
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